Predictions in eternal inflation 
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In generic models of cosmological inflation, quantum fluctuations strongly influence the spacetime 
metric and produce infinitely many regions where the end of inflation (reheating) is delayed until 
arbitrarily late times. The geometry of the resulting spacetime is highly inhomogeneous on scales 
of many Hubble sizes. The recently developed string-theoretic picture of the "landscape" presents a 
similar structure, where an infinite number of de Sitter, flat, and anti-de Sitter universes are nucle- 
ated via quantum tunneling. Since observers on the Earth have no information about their location 
within the eternally inflating universe, the main question in this context is to obtain statistical pre- 
dictions for quantities observed at a random location. I describe the problems arising within this 
statistical framework, such as the need for a volume cutoff and the dependence of cutoff schemes 
on time slicing and on the initial conditions. After reviewing different approaches and mathemat- 
ical techniques developed in the past two decades for studying these issues, I discuss the existing 
proposals for extracting predictions and give examples of their applications. 



I. ETERNAL INFLATION 

The general idea of eternally inflating spacetime was 
first introduced and developed in the 1980s [l|, 0, H, 0] m 
the context of slow-roll inflation. Let us begin by review- 
ing the main features of eternal inflation, following these 
early works. 

A prototypical model contains a minimally coupled 
scalar field <f> (the "inflaton") with an effective potential 
V(4>) that is sufficiently flat in some range of <f>. When 
the field <f> has values in this range, the spacetime is ap- 
proximately de Sitter with the Hubble rate 
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V{ct>) = H(<j>). 



(1) 



(We work in units where G = c = h = 1.) The value of H 
remains approximately constant on timescales of several 
Hubble times (At > H~ r ), while the field <j> follows the 
slow-roll trajectory <fi SI (t). Quantum fluctuations of the 
scalar field <t> in de Sitter background grow linearly with 
time 1,0,0], 
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at least for time intervals At of order several H _ 1 . Due to 
the quasi-exponential expansion of spacetime during in- 
flation, Fourier modes of the field <j) are quickly stretched 
to super-Hubble length scales. However, quantum fluc- 
tuations with super-Hubble wavelengths cannot main- 
tain Cjuantum coherence and become essentially classi- 
cal [1,0, B B EH ; this issue is discussed in more detail in 
Sec. II Bl below. The resulting field evolution <j>(t) can be 
visualized 0,H,[lH as a Brownian motion with a "random 
jump" of typical step size A(f> ~ H / (2tt) during a time 
interval At ~ H~ x , superimposed onto the determinis- 
tic slow-roll trajectory 4> SI (t). A statistical description 
of this "random walk"-type evolution <p(t) is reviewed in 

seeing 

The "jumps" at points separated in space by many 
Hubble distances are essentially uncorrelated; this is an- 
other manifestation of the well-known "no-hair" property 



of de Sitter space Thus the field 4> becomes 

extremely inhomogeneous on large (super-horizon) scales 
after many Hubble times. Moreover, in the semi-classical 
picture it is assumed that the local expansion rate 
a/a = H(4>) tracks the local value of the field (p(t,x) ac- 
cording to the Einstein equation (JTJ) . Here a(t, x) is the 
scale factor function which varies with x only on super- 
Hubble scales, a(t,x)Ax > H~ l . Hence, the spacetime 
metric can be visualized as having a slowly varying, "lo- 
cally de Sitter" form (with spatially flat coordinates x) , 



dx^da 
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The deterministic trajectory cj> SI (t) eventually reaches 
a (model-dependent) value 0* signifying the end of the 
slow-roll inflationary regime and the beginning of the re- 
heating epoch (thermalization). Since the random walk 
process will lead the value of </> away from cj> = in 
some regions, reheating will not begin everywhere at 
the same time. Moreover, regions where <j> remains in 
the inflationary range will typically expand faster than 
regions near the end of inflation where V(4>) becomes 
small. Therefore, a delay of the onset of reheating will 
be rewarded by additional expansion of the proper 3- 
volume, thus generating more regions that are still in- 
flating. This feature is called "self-reproduction" of the 
inflationary spacetime Q . Since each Hubble-size region 
evolves independently of other such regions, one may vi- 
sualize the spacetime as an ensemble of inflating Hubble- 
size domains (Fig. [T]) . 

The process of self-reproduction will never result in 
a global reheating if the probability of jumping away 
from (f> = 0» and the corresponding additional volume 
expansion factors are sufficiently large. The correspond- 
ing quantitative conditions and their realization in typical 
models of inflation are reviewed in Sec. lIII Al Under these 
conditions, the process of self-reproduction of inflating re- 
gions continues forever. At the same time, every given 
comoving worldlinc (except for a set of measure zero; see 
Sec. IIII A|) will sooner or later reach the value <fi = 0* 
and enter the reheating epoch. The resulting situation is 
known as "eternal inflation" 3] . More precisely, the term 
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Figure 1: A qualitative diagram of self-reproduction during 
inflation. Shaded spacelike domains represent Hubble-size re- 
gions with different values of the inflaton field cf>. The time 
step is of order H -1 . Dark-colored shades are regions under- 
going reheating (cS = </>*); lighter-colored shades are regions 
where inflation continues. On average, the number of inflating 
regions grows with time. 




Figure 2: A 1+1-dimensional slice of the spacetime struc- 
ture in an eternally inflating universe (numerical simulation 
in Ref. [3). Shades of different color represent different, 
causally disconnected regions where reheating took place. The 
reheating surface is the line separating the white (inflating) 
domain and the shaded domains. 



"eternal inflation" means future-eternal self-reproduction 
of inflating regions [la]. To emphasize the fact that self- 
reproduction is due to random fluctuations of a field, one 
refers to this scenario as "eternal inflation of random-walk 
type." Below we use the terms "eternal self-reproduction" 
and "eternal inflation" interchangeably. 

Observers like us may appear only in regions where 
reheating already took place. Hence, it is useful to con- 
sider the locus of all reheating events in the entire space- 
time; in the presently considered example, it is the set 
of spacetime points x there 4>{x) — 0*. This locus is 
called the reheating surface and is a noncompact, space- 
like three-dimensional hypersurface (lfjl . [l8T ] . It is impor- 
tant to realize that a finite, initially inflating 3-volume 
of space may give rise to a reheating surface having an 
infinite 3-volume, and even to infinitely many causally 



disconnected pieces of the reheating surface, each hav- 
ing an infinite 3-volume (see Fig. [2] for an illustration). 
This feature of eternal inflation is at the root of several 
technical and conceptual difficulties, as will be discussed 
below. 

Everywhere along the reheating surface, the reheating 
process is expected to provide appropriate initial condi- 
tions for the standard "hot big bang" cosmological evolu- 
tion, including nucleosynthesis and structure formation. 
In other words, the reheating surface may be visualized 
as the locus of the "hot big bang" events in the space- 
time. It is thus natural to view the reheating surface as 
the initial equal-time surface for astrophysical observa- 
tions in the post-inflationary epoch. Note that the ob- 
servationally relevant range of the primordial spectrum 
of density fluctuations is generated only during the last 
60 e-foldings of inflation. Hence, the duration of the in- 
flationary epoch that preceded reheating is not directly 
measurable beyond the last 60 e-foldings; the total num- 
ber of e-foldings can vary along the reheating surface and 
can be in principle arbitrarily large. 2 

The phenomenon of eternal inflation is also found in 
multi-field models of inflation [22|, l23l. as well as in sce- 
narios based on Brans-Dicke theory [23, EE 26] , topologi- 
cal inflation 2?], .28;], braneworld inflation [23], "recycling 
universe" [30], and the string theory landscape 31j. In 
some of these models, quantum tunneling processes may 
generate "bubbles" of a different phase of the vacuum (see 
Sec. Ill El for more details). Bubbles will be created ran- 
domly at various places and times, with a fixed rate per 
unit 4- volume. In the interior of some bubbles, additional 
inflation may take place, followed by a new reheating sur- 
face. The interior structure of such bubbles is sketched in 
Fig. [3J The nucleation event and the formation of bub- 
ble walls is followed by a period of additional inflation, 
which terminates by reheating. Standard cosmological 
evolution and structure formation eventually give way to 
a A-dominated universe. Infinitely many galaxies and 
possible civilizations may appear within a thin spacelike 
slab running along the interior reheating surface. This 
reheating surface appears to interior observers as an in- 
finite, spacelike hypersurface [13]. For this reason, such 
bubbles are called "pocket universes," while the spacetime 
is called a "multiverse." (Generally, the term "pocket uni- 
verse" refers to a noncompact, connected component of 
the reheating surface [33|.) 

In scenarios of this type, each bubble is causally dis- 
connected from most other bubbles. 3 Hence, bubble nu- 



It is worth emphasizing that the term "eternal inflation" refers 
to future-eternity of inflation in the sense described above, but 
does not imply past-eternity. In fact, inflationary spacetimes are 
generically not past-eternal Il6l. ll7l| . 



For instance, it was shown that holographic considerations do 
not place any bounds on the total number of e-foldings during 
inflation |2(l . For recent attempts to limit the number of e- 
foldings using a different approach, see e.g. [2l| . Note also that 
the effects of "random jumps" are negligible during the last 60 
e-foldings of inflation, since the produced perturbations must be 
of order 10 -5 according to observations. 

Collisions between bubbles are rare 1341: however, effects of bub- 
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Figure 3: A spacetime diagram of a bubble interior. The in- 
finite, spacelike reheating surface is shown in darker shade. 
Galaxy formation is possible within the spacetime region in- 
dicated. 

cleation events may generate infinitely many statistically 
inequivalcnt, causally disconnected patches of the reheat- 
ing surface, every patch giving rise to a possibly infinite 
number of galaxies and observers. This feature signifi- 
cantly complicates the task of extracting physical predic- 
tions from these models. This class of models is referred 
to as "eternal inflation of tunneling type." 

In the following subsections, I discuss the motivation 
for studying eternal inflation as well as physical justifica- 
tions for adopting the effective stochastic picture. Differ- 
ent techniques developed for describing eternal inflation 
are reviewed in Sec. [TTJ Section IIIII contains an overview 
of methods for extracting predictions and a discussion of 
the accompanying "measure problem." 



A. Some motivation 

The hypothesis of cosmological inflation was invoked 
to explain several outstanding puzzles in observational 
data [36j | . However, some observed quantities (such as the 
cosmological constant A or elementary particle masses) 
may be expectation values of slowly-varying effective 
fields Xa- Within the phenomcnological approach, we are 
compelled to consider also the fluctuations of the fields 
Xa during inflation, on the same footing as the fluctua- 
tions of the inflaton 4>. Hence, in a generic scenario of 
eternal inflation, all the fields Xa arrive at the reheating 
surface <f> — (/)» with values that can be determined only 
statistically. Observers appearing at different points in 
space may thus measure different values of the cosmo- 
logical constant, elementary particle masses, spectra of 
primordial density fluctuations, and other cosmological 
parameters. 

It is important to note that inhomogeneities in observ- 
able quantities are created on scales far exceeding the 
Hubble horizon scale. Such inhomogeneities are not di- 



ble collisions are observable in principle |35t . 



rectly accessible to astrophysical experiments. Neverthe- 
less, the study of the global structure of eternally inflat- 
ing spacetime is not merely of academic interest. Funda- 
mental questions regarding the cosmological singularities, 
the beginning of the Universe and of its ultimate fate, as 
well as the issue of the cosmological initial conditions 
all depend on knowledge of the global structure of the 
spacetime as predicted by the theory, whether or not this 
global structure is directly observable (see e.g. [37l l38l|). 
In other words, the fact that some theories predict eter- 
nal inflation influences our assessment of the viability of 
these theories. In particular, the problem of initial con- 
ditions for inflation [3^ | is significantly alleviated when 
eternal inflation is present. For instance, it was noted 
early on that the presence of eternal self-reproduction 
in the "chaotic" inflationary scenario [Io| essentially re- 
moves the need for the fine-tuning of the initial condi- 
tions More recently, constraints on initial condi- 
tions were studied in the context of self-reproduction in 
models of quintessence [42[ and fc-inflation [43| . 

Since the values of the observable parameters Xa are 
random, it is natural to ask for the probability distribu- 
tion of Xa that would be measured by a randomly chosen 
observer. Understandably, this question has been the 
main theme of much of the work on eternal inflation. 
Obtaining an answer to this question promises to estab- 
lish a more direct contact between scenarios of eternal 
inflation and experiment. For instance, if the probability 
distribution for the cosmological constant A were peaked 
near the experimentally observed, puzzlingly small value 
(see e.g. [44J for a review of the cosmological constant 
problem) , the smallness of A would be explained as due 
to observer selection effects rather than to fundamen- 
tal physics. Considerations of this sort necessarily in- 
volve some anthropic reasoning; however, the relevant 
assumptions are minimal. The basic goal of theoretical 
cosmology is to select physical theories of the early uni- 
verse that are most compatible with astrophysical obser- 
vations, including the observation of our existence. It ap- 
pears reasonable to assume that the civilization of Planet 
Earth evolved near a randomly chosen star compatible 
with the development of life, within a randomly chosen 
galaxy where such stars exist. Many models of inflation 
generically include eternal inflation and hence predict the 
formation of infinitely many galaxies where civilizations 
like ours may develop. It is then also reasonable to as- 
sume that our civilization is typical among all the civi- 
lizations that evolved in galaxies formed at any time in 
the universe. This assumption is called the "principle of 
mediocrity" [181 ]. 

To use the "principle of mediocrity" for extracting sta- 
tistical predictions from a model of eternal inflation, one 
proceeds as follows [HI, EE} . In the example with the 
fields Xa described above, the question is to determine 
the probability distribution for the values of Xa that a 
random observer will measure. Presumably, the values 
of the fields Xa do not directly influence the emergence 
of intelligent life on planets, although they may affect 
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the efficiency of structure formation or nucleosynthesis. 
Therefore, we may assume a fixed, x a -dependent mean 
number of civilizations f c iv(Xa) per galaxy and proceed 
to ask for the probability distribution Pg(Xo) of Xa near 
a randomly chosen galaxy. The observed probability dis- 
tribution of Xa will then be 

P(Xa) = PoiXaWiAXa)- (4) 

One may use the standard "hot big bang" cosmology to 
determine the average number va(Xa) of suitable galax- 
ies per unit volume in a region where reheating occurred 
with given values of Xa', in any case, this task does not 
appear to pose difficulties of principle. Then the compu- 
tation of Pg{Xo) is reduced to determining the volume- 
weighted probability distribution V(xa) for the fields Xa 
within a randomly chosen 3-volume along the reheating 
surface. The probability distribution of Xa will be ex- 
pressed as 

P{Xa) = V(Xa)v G (Xa)v c i v (Xa)- (5) 

However, defining V(xa) turns out to be far from straight- 
forward since the reheating surface in eternal inflation is 
an infinite 3-surface with a complicated geometry and 
topology. The lack of a natural, unambiguous, unbi- 
ased measure on the infinite reheating surface is known 
as the "measure problem" in eternal inflation. Existing 
approaches and measure prescriptions are discussed in 
Sec. IIIIl where two main alternatives (the "volume-based" 
and "worldline-based" measures) are presented. In Sec- 
tions IIIIBI and IIIIDI I give arguments in favor of using 
the volume-based measure for computing the probabil- 
ity distribution of values Xa measured by a random ob- 
server. The volume-based measure has been applied to 
obtain statistical predictions for the gravitational con- 
stant in Brans-Dicke theories I24T. [2511 . cosmological con- 
stant (dark energy) |46l l47l l48l . l49l . [HoL HH, particle 
physics parameters |52l. l53l. l54j. and the amplitude of 
primordial density perturbations [H, [5l|, [H, [56| . 

The issue of statistical predictions has recently come to 
the fore in conjunction with the discovery of the string 
theory landscape. According to various estimates, one 
expects to have between 10 500 and 10 1500 possible vac- 
uum states of string theory [3l], H3, [H, HjJ The 
string vacua differ in the geometry of spacetime compact- 
ification and have different values of the effective cosmo- 
logical constant (or "dark energy" density). Transitions 
between vacua may happen via the well-known Coleman- 
deLuccia tunneling mechanism [32| . Once the dark en- 
ergy dominates in a given region, the spacetime becomes 
locally de Sitter. Then the tunneling process will cre- 
ate infinitely many disconnected "daughter" bubbles of 
other vacua. Observers like us may appear within any 
of the habitable bubbles. Since the fundamental theory 
does not specify a single "preferred" vacuum, it remains 
to try determining the probability distribution of vacua 
as found by a randomly chosen observer. The "volume- 
based" and "worldline-based" measures can be extended 



to scenarios with multiple bubbles, as discussed in more 
detail in Sec. IIIIDI Some recent results obtained using 
these measures are reported in Refs. [gj Ir32l [63|. 



B. Physical justifications of the semiclassical 
picture 

The standard framework of inflationary cosmology as- 
serts that vacuum quantum fluctuations with super- 
horizon wavelengths become classical inhomogeneities of 
the field <p. The calculations of cosmological density per- 
turbations generated during inflation 0, [HI, [U, [65|, [6|| 
H3, HI] also assume that a "classicalization" of quantum 
fluctuations takes place via the same mechanism. In the 
calculations, the statistical average {8(f) 2 ) of classical fluc- 
tuations on super-Hubble scales is simply set equal to the 
quantum expectation value (0| <fi 2 |0) in a suitable vac- 
uum state. While this approach is widely accepted in 
the cosmology literature, a growing body of research is 
devoted to the analysis of the quantum-to-classical tran- 
sition during inflation (see e.g. (69j for an early review). 
Since a detailed analysis would be beyond the scope of 
the present text, I merely outline the main ideas and ar- 
guments relevant to this issue. 

A standard phenomenological explanation of the "clas- 
sicalization" of the perturbations is as follows. For sim- 
plicity, let us restrict our attention to a slow-roll infla- 
tionary scenario with one scalar field <f>. In the slow roll 
regime, one can approximately regard as a massless 
scalar field in de Sitter background spacetime [1| . Due to 
the exponentially fast expansion of de Sitter spacetime, 
super-horizon Fourier modes of the field (f> are in squeezed 
quantum states with ex p onentiall y la rge (~ e Ht ) squeez- 
ing parameters [zO, [zE [zl M, Izl tz| • Such highly 
squeezed states have a macroscopically large uncertainty 
in the field value <j> and thus quickly decohere due to 
interactions with gravity and with other fields. The re- 
sulting mixed state is effectively equivalent to a statis- 
tical ensemble with a Gaussian distributed value of <j>. 
Therefore one may compute the statistical average (<5</> 2 ) 

as the quantum expectation value (0| (f> 2 |0) and interpret 
the fluctuation 8<f> as a classical "noise." A heuristic de- 
scription of the "classicalization" [i[ is that the quantum 
commutators of the creation and annihilation operators 
of the field modes, [a, a)\ = 1, are much smaller than the 
expectation values (a^a) 3> 1 and are thus negligible. 

A related issue is the backreaction of fluctuations of 
the scalar field <fi on the metric. 4 According to the stan- 
dard theory (see e.g. [H, [H[ for reviews), the perturba- 
tions of the metric arising due to fluctuations of (f> are 
described by an auxiliary scalar field (sometimes called 



4 The backreaction effects of the long-wavelength fluctuations of 
a scalar field du ring inflation have been investigated extensively 
(see e.g. m H H EL Ek EE HI HI)- 
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the "Sasaki-Mukhanov variable") in a fixed de Sitter back- 
ground. Thus, the "classicalization" effect should apply 
equally to the fluctuations of <fr and to the induced metric 
perturbations. At the same time, these metric perturba- 
tions can be viewed, in an appropriate coordinate system, 
as fluctuations of the local expansion rate H(4>) due to 
local fluctuations of <f> [3, 0, IM] • Thus one arrives at the 
picture of a "locally de Sitter" spacetime with the met- 
ric ([3]), where the Hubble rate a/a = H{cf)) fluctuates on 
super-horizon length scales and locally follows the value 
of <fi via the classical Einstein equation ([lj . 

The picture as outlined is phenomenological and does 
not provide a description of the quantum-to-classical 
transition in the metric perturbations at the level of field 
theory. For instance, a fluctuation of <j> leading to a lo- 
cal increase of H(S) necessarily violates the null energy 
condition [M H3, [M] . The cosmological implications of 
such "semiclassical" fluctuations (see e.g. the scenario of 
"island cosmology" [89|, |9(J, [9l|) cannot be understood 
in detail within the framework of the phenomenological 
picture. 

A more fundamental approach to describing the 
quantum-to-classical transition of perturbations was de- 
veloped using non-equilibrium quantum field theory and 
the influence functional formalism [M IM IM IM] ■ In this 
approach, decoherence of a pure quantum state of <f> into 
a mixed state is entirely due to the self-interaction of the 
field (j>. In particular, it is predicted that no decoherence 
would occur for a free field with V(<p) — \m§ 2 . This 
result is at variance with the accepted paradigm of "clas- 
sicalization" as outlined above. If the source of the "noise" 
is the coupling between different perturbation modes of <j>, 
the typical amplitude of the "noise" will be second-order 
in the perturbation. This is several orders of magnitude 
smaller than the amplitude of "noise" found in the stan- 
dard approach. Accordingly, it is claimed [M [93] that 
the magnitude of cosmological perturbations generated 
by inflation is several orders of magnitude smaller than 
the results currently accepted as standard, and that the 
shape of the perturbation spectrum depends on the de- 
tails of the process of "classicalization" [98(. Thus, the 
results obtained via the influence functional techniques 
do not appear to reproduce the phenomenological pic- 
ture of "classicalization" as outlined above. This mis- 
match emphasizes the need for a deeper understanding 
of the nature of the quantum-to-classical transition for 
cosmological perturbations. 

Finally, let us mention a different line of work which 
supports the "clas sicalization" picture. In Refs. [1, l99l . 
llOOl . llOll Il02l . Il03| ] . calculations of (renormalized) expec- 
tation values such as (<^ 2 ), (0 4 ), etc., were performed for 
field operators in a fixed de Sitter background. The 
results were compared with the statistical averages 

J P{(f>,t)<t> 2 d<f>, J P((f),t)(t) 4 d(j), etc., (6) 



Sec. Ill Ap . It was shown that the leading late-time 
asymptotics of the quantum expectation values coincide 
with the corresponding statistical averages ©. These 
results appear to validate the "random walk" approach, 
albeit in a limited context (in the absence of backreac- 
tion) . 



II. STOCHASTIC APPROACH TO INFLATION 

The stochastic approach to inflation is a semiclassi- 
cal, statistical description of the spacetime resulting from 
quantum fluctuations of the in flaton field(s) and their 
backreaction on the metric [J, i, EM EM EM, EM EM 
EM IllCt Illll Ill2j |. In this description, the spacetime 
remains everywhere classical but its geometry is deter- 
mined by a stochastic process. In the next subsections I 
review the main tools used in the stochastic approach for 
calculations in the context of random- walk type, slow-roll 
inflation. Models involving tunncling-typc eternal infla- 
tion are considered in Sec. Ill El 



A. Random walk-type eternal inflation 

The important features of random walk-type eternal 
inflation can be understood by considering a simple slow- 
roll inflationary model with a single scalar field <f> and a 
potential V(<fi). The slow-roll evolution equation is 

• 1 dV 1 dH rj _. _ 

where H((f>) is defined by Eq. JT]) and v((f>) is a model- 
dependent function describing the "velocity" </> of the de- 
terministic evolution of the field <fi. The slow-roll trajec- 
tor y </>sr ( iQ, w hich is a solution of Eq. ([7]), is an attrac- 
tor [1131 . Ill4j for trajectories starting with a wide range 
of initial conditions. 5 

As discussed in Sec. II Bl the super-horizon modes of 
the field <j) are assumed to undergo a rapid quantum- 
to-classical transition. Therefore one regards the spatial 
average of <fi on scales of several H~ x as a classical field 
variable. The spatial averaging can be described with 
help of a suitable window function, 

(0(x)} = J W(x - y)0(y)d 3 y- (8) 

It is implied that the window function W(x) decays 
quickly on physical distances a |x| of order several H~ l . 
From now on, let us denote the volume-averaged field 
simply by <f> (no other field <j) will be used). 



where the distribution P{4>, t) describes the "random 
walk" of the field <f> in the Fokker-Plank approach (see 



5 See Ref. [43l | for a precise definition of an attractor trajectory in 
the context of inflation. 
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As discussed above, the influence of quantum fluctu- 
ations leads to random "jumps" superimposed on top of 
the deterministic evolution of the volume-averaged field 
</>(i,x). This may be described by a Langevin equation 
of the form 0] 



0(i,x) = v(</>)+N(t,x), 



(9) 



where N(t,x) stands for "noise" and is assumed to be a 
Gaussian random function with known correlator 0, Ill5l 

EIlQTi} 



(N(t,x)N(t,k)) = C(t,t, 



(10) 



An explicit form of the correlator C depends on the spe- 
cific window funct ion W used for averaging the field <fi 
on Hubble scales [ll7j . However, the window function 
W is merely a phenomenological device used in lieu of a 
complete ab initio derivation of the stochastic inflation 
picture. One expects, therefore, that results of calcula- 
tions should be robust with respect to the choice of W. 
In other words, any uncertainty due to the choice of the 
window function must be regarded as an imprecision in- 
herent in the method. For instance, a robust result in 
this sense is an exponentially fast decay of correlations 
on time scales At > H^ 1 , 

C(t,f,|x-x|;^) acexp(-2H(4>)\t-t\) , (11) 

which holds for a wide class of window functions [1171 ]. 

For the purposes of the present consideration, we only 
need to track the evolution of <f>(t, x) along a single co- 
moving worldline x = const. Thus, we will not need an 
explicit form of C(t, t, |x — x| ; <f>) but merely the value at 
coincident points t = t, x = x, which is computed in the 
slow-roll inflationary scenario as Q 



C( t> t,0;<P) = ^- 



(12) 



(This represents the fluctuation ^ accumulated during 
one Hubble time, At = H^ 1 .) Due to the property (fTTj) . 
one may neglect correlations on time scales At > H^ 1 in 
the "noise" field. 6 Thus, the evolution of <f> on time scales 
At > H^ 1 can be described by a finite-difference form of 
the Langevin equation ©, 



<f>(t + At) - 4>{t) = v(<f>)At + v / 2D(c/>)At^(t), (13) 
where 

H 3 (cj>) 



D{4>) 



(14) 



Taking these cor relations into account leads to a picture of "color 
noise" |l!8l . Ilisjl . In what follows, we only consider the simpler 
picture of "white noise" as an approximation adequate for the 
issues at hand. 



and £ is a normalized random variable representing 
"white noise," 



(0=0, (£ 2 ) = 1, 
+ At)) = for At > H~ 



(15) 
(16) 



Equation ()13|) is interpreted as describing a Brownian 
motion cf>(t) with the systematic "drift" v(<f>) and the "dif- 
fusion coefficient" D(<fi). In a typical slow-roll inflation- 
ary scenario, there will be a range of <fi where the noise 
dominates over the deterministic drift, 



v(cj))At < y/2D(<j))At, At = H~ 



(17) 



Such a range of 4> is called the "diffusion-dominated 
regime." For <j> near the end of inflation, the amplitude 
of the noise is very small, and so the opposite inequality 
holds. This is the "deterministic regime" where the ran- 
dom jumps can be neglected and the field <fi follows the 
slow-roll trajectory. 



B. Fokker-Planck equations 

A useful description of the statistical properties of 
4>(t) is furnished by the probability density P(4>,t)d<f) 
of having a value <j> at time t. As in the case of the 
Langevin equation, the values (f>(t) are measured along a 
single, randomly chosen comoving worldline x = const. 
The probability distribution P(<p, t) satisfies the Fokker- 
Planc k (FP) equation whose standard derivation we 

omit [imfm] . 



d t P = [~v(cp)P + d {D(<t>)P)] 



(18) 



The coefficients v((j)) and D{<p) are in general model- 
dependent and need to be calculated in each particular 
scenario. These calculations require only the knowledge 
of the slow-roll trajectory and the mode functions of the 
quantized scalar perturbations. For ordinary slow-roll 
inflation with an effective potential V(<j>), the results are 
well-known expressions (O and (fT4]) . The correspond- 
ing expressions for models of fc-inflation were derived in 
Ref. 14311 u sing the relevant quantum theory of perturba- 
tions jl22l |. 

It is well known that there exists a "factor ordering" 
ambiguity in translating the Langevin equation into the 
FP equation if the amplitude of the "noise" depends on 
the position. Specifically, the factor D{(j>) in Eq. (fl3)) 
may be replaced by D(<f> + 6 At), where < 8 < 1 is an 
arbitrary constant. With 9^0, the term (DP) in 
Eq. (TT5]) will be replaced by a different ordering of the 
factors, 



d H (DP) -f d (D^P)] 



(19) 



Popular choices 8 — and 8 — | are called the Ito and 
the Stratonovich factor ordering respectively. Motivated 
by the considerations of Ref. [l23j | , we choose 8 — as 
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shown in Eqs. (fT5| and (fT5|) . Given the phenomenolog- 
ical nature of the Langevin equation (|13[) . one expects 
that any ambiguity due to the choice of 9 represents an 
imprecision inherent in the stochastic a ppro ach. This 
imprecision is typically of order H 2 <C 1 |l24j . 

The quantity P(<p, t) may be also interpreted as the 
fraction of the comoving volume (i.e. coordinate volume 
c? 3 x) occupied by the field value (f> at time t. Another 
important characteristic is the volume-weighted distribu- 
tion Pv{4>,t)d(j) 1 which is defined as the proper 3-volumc 
(as opposed to the comoving volume) of regions having 
the value <fr at time t. (To avoid considering infinite vol- 
umes, one may restrict one's attention to a finite comov- 
ing domain in the universe and normalize Pv(<P, t) to unit 
volume at some initial time t = to.) The volume d istri- 
bution satisfies a modified FP equation 0, llOTt Il09j ] , 

d t Pv = [-v(cj>)P v + c> (D(cj))P v )] + 3H(4>)P V , (20) 

which differs from Eq. (fl~8|) by the term 3HPy that de- 
scribes the exponential growth of 3-volume in inflating 
regions. 7 

Presently we consider scenarios with a single scalar 
field; however, the formalism of FP equations can be 
straightforwardly extended to multi-field models (see 
e.g. Ref. |l23l |). For instance, the FP equation for a two- 
field model is 

d t p = a 00 (DP)+d xx (DP)-d^ {v^p)-d x (v x p) , (21) 

where D(<f>, %), v,p((j), x)j and v x ((j>, x) are appropriate co- 
efficients. 

C. Methods of solution 

In principle, one can solve the FP equations forward in 
time by a numerical method, starting from a given initial 
distribution at t = to- To specify the solution uniquely, 
the FP equations must be supplemented by boundary 
cond itions at both ends of the inflating range of <f> [llll 
Ill2| . At the reheating boundary (</> = 0*), one imposes 
the "exit" boundary conditions, 

fy[^WW=o, fyTO)JV]^. = o. (22) 

These boundary conditions express the fact that random 
jumps are very small at the end of inflation and cannot 
move the value of </> away from <j> = <^>*. If the potential 
V((j)) reaches Planck energy scales at some <j> = ^max (this 
happens generally in "chaotic" type inflationary scenarios 
with unbounded potentials), the semiclassical picture of 
spacetime breaks down for regions with (f> ~ </> m ax- Hence, 



7 A more formal derivation of Eq. (1 20 I t as well as details of the in- 
terpretation of the distributions P and Py in terms of ensembles 
of worldlines can be found in Ref. [125| . 



a boundary condition must be imposed also at (f> = </> max - 
For instance, one can use the absorbing boundary condi- 
tion, 

P(0max) = 0, (23) 

which means that Planck-ene rgy region s with <f> — max 
disappear from consideration [llll Ill2j . 

Once the boundary conditions are specified, one may 
write the general solution of the FP equation (fT5|) as 

P(4>,t) = Y / CxP W ^)e Xt , (24) 
x 

where the sum is performed over all the eigenvalues A of 
the differential operator 

LP = d^[~v(cj))P + d^(D(cj))P)}, (25) 

and the corresponding eigenfunctions P^ are defined by 

LP^ (0) = AP( A) (0). (26) 

The constants C\ can be expressed through the initial 
distribution P(4>,to). 

By an appropriate change of variables <fi — > z, P(4>) — > 
F(z), the operat or L may be brought into a manifestly 
self-adjoint form pM [Tol 0M QM EH, EH , 

L^—z+U(z). (27) 

Then one can show that all the eigenvalues A of L are non- 
positive; in particular, the (algebraically) largest eigen- 
value A max = —7 < is nondegenerate and the cor- 
respondi ng eig enfunction p( Amax )(0) is everywhere pos- 
itive [43|, Il24j . Hence, this eigenfunction describes the 
late-time asymptotic of the distribution P{(j) 1 t) 1 

PliiJaF^-'We- 1 ' (28) 

The distribution P^™ 1 ^' ((f) is the "stationary" distribu- 
tion of <f> per comoving volume at late times. The expo- 
nential decay of the distribution P(<p,t) means that at 
late times most of the comoving volume (except for an 
exponentially small fraction) has finished inflation and 
entered reheating. 

Similarly, one can represent the general solution of 
Eq. j20]) by 

Pv(^,t) = ^C^ A) (0)e At , (29) 

A 

where 

[L + 3#(0)]P( A) (0) = AP( A )(0). (30) 

By the same method as for the operator L, it is possible 
to show that the spectrum of eigenvalues A of the op- 
erator L + 3H(<fi) is bounded from above and that the 
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largest eigenvalue A max = 7 admits a nondegenerate, ev- 
erywhere positive eigenfunction P^\(f). However, the 
largest eigenvalue 7 may be either positive or negative. 
If 7 > 0, the late-time behavior of Pv(4>, t) is 



(31) 



which means that the total proper volume of all the in- 
flating regions grows with time. This is the behavior ex- 
pected in eternal inflation: the number of independently 
inflating domains increases without limit. Thus, the con- 
dition 7 > is the criterion for the presence of eternal 
self-reproduction of inflating domains. The correspond- 
ing distrib ution P^((j>) is called the "stationary" distri- 
bution [H, EE El Ezl • 

If 7 < 0, eternal inflation does not occur and the entire 
space almost surely (i.e. with probability 1) enters the 
reheating epoch at a finite time. 

If the poten t ial V {(f) is of "new" inflationary type @, 

1 1281 . I129L Il30l Il3l| and has a global maximum at say 
(f> = (j>o, the eigenvalues 7 and 7 can be es tima ted (under 
the usual slow-roll assumptions on V) as [l24j 

^~ r'vfl\ H ^ <0 ' 7«3ff(^o)>0. (32) 

Therefore, eternal inflation is generic in the "new" infla- 
tionary scenario. 

Let us comment on the possibility of obtaining solu- 
tions P((f),t) in practice. With the potential V{<f>) — 
A</> 4 , the full time-dependent FP equation (jT5J) can be 
solved ana lytically via a nonlinear change of variable 
^ -> (j)~ 2 [lH, [[Hi [l33[. This exact solution, as well 
as an approximate solution P((f>, t) for a general poten- 
tial, can be also obtained using the saddle-point evalu- 
ation of a path- integral expression for P(<f>,t) [l34| . In 
some cases the eigenvalue equation LP^ = \P^ may 
be reduced to an exactly solvable Schrodinger equation. 
These cases include potentials of the form V((f>) = A e^, 



V(<p) = \(f>~' 2 , V(4>) = Acosh^(/x</>); see e.g. [L2j for 
other examples. 

A general approximat e me t hod for d etermining P(</>, t) 
for arbitrary potentials [l35l Il36l . Il37l | consists of a per- 
turbative expansion, 



ct>(t) = Mt) + $Mt) + 6Mt) 



(33) 



applied directly to the Langevin equation. The result is 
(at the lowest order) a Gaussian approximation with a 
time-dependent mean and variance [135j |. 



P(4>,t)* 

a 2 (t)^ 



1 



H' 2 (<j> SI 

7T 



exp 



2a 2 (t) 



#73 ' 



b (t)=cp sl .(t) + ^a 2 (t) + ^- 



4tt 



H 3 

111 

W 2 



H' 2 



(34) 

(35) 
(36) 



where <j> SI (i) is the slow-roll trajectory and is the ini- 
tial value of </>. While methods based on the Langevin 
equation do not take into account boundary conditions 
or volume weighting effects, the formula (|34[) provides an 
adequate approximation to the distribution P(<p,t) in a 
useful range of <j> and t [l37j | . 



D. Gauge dependence issues 

An important feature of the FP equations is their de- 
pendence on the choice of the time variable. One can 
consider a replacement of the form 



t 



t, dr — T(<fi)dt, 



(37) 



understood in the sense of integrating along comoving 
worldlines x = const, where T(<p) > is an arbitrary 
function of the field. For instance, a possible choice is 
T((fi) = H(<f>), which makes the new time variable di- 
mcnsionless, 



Hdt = In a. 



(38) 



This time variable is called "scale factor time" or "e- 
folding time" since it measures the number of e-foldings 
along a comoving worldlinc. 

The distributions P(<f>, r) and Pv(<P,t) are defined as 
before, except for considering the 3-volumes along hy- 
persurfaces of equal r. These distributions satisfy FP 
equations similar to Eqs. P^)) - PD|) . With the replace- 
ment (|37p. the coe fficien ts of the new FP equations are 
modified as follows jl24| . 



D(0) 



D(<f>) 



v{<t>) 



y(4>) 
T(4,y 



(39) 



while the "growth" term 3HPy in Eq. (f20|) is replaced 
by 3HT~ 1 Pv ■ The change in the coefficients may sig- 
nificantly alter the qualitative behavior of the solutions 
of the FP equations. For instance, stationary distribu- 
tions defined through the proper time t and the e-folding 
time r = I n a were found to have radically different be- 
havior [l8l [TT^ . [l27j . This sensitivity to the choice of the 
"time gauge" r is unavoidable since hypersurfaces of equal 
r may preferentially select regions with certain proper- 
ties. For instance, most of the proper volume in equal- 
t hypersurfaces is filled with regions that have gained 
expansion by remaining near the top of the potential 
V{4>), while hypersurfaces of equal scale factor will under- 
represent those regions. Thus, a statement such as "most 
of the volume in the Universe has values of <j> with high 
V{4>)" is largely gauge-dependent. 

In the early works on eternal inflation (2f| lllll . El 

Il27l |. the late-time asymptotic distribution of volume 
Py\<j)) along hypersurfaces of equal proper time [see 
Eq. (|3ip ] was interpreted as the stationary distribution 
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of field values in the universe. However, the high sensi- 
tivity of this distribution to the choice of the time vari- 
able make s thi s interpretation unsatisfactory. Also, it 
was noted [l38j | that equal-proper time volume distribu- 
tions predict an unacceptably small probability for the 
currently observed CMB temperature. The reason for 
this result is the extreme bias of the proper-time gauge to- 
wards over-re pre s ent ing regions where reheating occurred 
very recently jl8l . Il39| . One might ask whether hypersur- 
faces of equal scale factor or some other choice of time 
gauge wo uld p rovide less biased answers. However, it 
turns out [125j that there exists no a priori choice of the 
time gauge r that provides unbiased equal-r probability 
distributions for all potentials V((f>) in models of slow-roll 
inflation (see Sec. IHICl for details). 

Although the FP equations necessarily involve a de- 
pendence on gauge, they do provide a useful statisti- 
cal picture of the distribution of fields in the universe. 
The FP techniques can also be used for deriving sev- 
eral gauge-independent results. For instance, the pres- 
ence of eternal inflation is a gauge- independent statement 
(see also Sec. IIII Aj) : if the largest eigenvalue 7 is posi- 
tive in one gaug e of the form (|3"T)) . then 7 > in ev- 
ery other gauge |l40| |. Using the FP approach, one can 
also c ompute th e fractal dimension of the inflating do- 
main [Ho, Il4l| and the probability of exiting inflation 
through a particular point 4>* of the reheating boundary 
in the configuration space (in case there exists more than 
one such point). 

The exit probability can be determined as follows [43l . 

1 1 23f | . Let us assume for simplicity that there are two 
possible exit points 0* and 4>e, and that the initial dis- 
tribution is concentrated at 6 — <t>n, i.e. 



boundary conditions 



p{4>,t = o) = 6{4>-<h), 



(40) 



where 4>e < 4>o < 4>*- The probability of exiting inflation 
through <f> = <f>E during a time interval [t, t + dt] is 



dp, 



-v(</> E )P(<l>E,t)dt 



(41) 



(note that v(4>e) < 0). Hence, the total probability of 
exiting through <j> = (j>E at any time is 



Pexit(0B) = / dp c ^ it (4> E ) = -v{4>e) I P((f)B,t)dt, 
Jo Jo 

(42) 

Introducing an auxiliary function F(<p) as 



F(</>) = -v((/>) P{cj>,t)dt 



(43) 



one can show that F(cj>) satisfies the gauge-invariant 
equation, 



,D 



= 0, df(-F) 
v 



0. 



(45) 



can be straightforwardly integrated and yields explicit 
expressions for the exit probability p e xit(</>_E) as a function 
of the initial value 4>o [43j . The exit probability p e xit (</>*) 
can be determined similarly. 



E. Self-reproduction of tunneling type 

Until now, we considered eternal self-reproduction due 
to random walk of a scalar field. Another important class 
of models includes self-reproduction due to bubble nucle- 
ation. 8 Suc h scenarios of eternal in flation were studied 
in Refs. P. [Mifl4l . [lli . [l4lL[M6t . 

In a locally de Sitter universe dominated by dark en- 
ergy, nucleation of bub bles of fal se vacuum may occur- 
due to tunneling 0, [H, EI3, EH ■ Since the bubble nu- 
cleation rate n per unit 4- volume is very small [H, 1 1491 ] , 



k = 0(l)H- 4 exp 



~ Sl IP. 



(46) 



where Si is the instanton action and H is the Hubble 
constant of the de Sitter background, bubbles will gener- 
ically not merge into a single false- vacuum domain [34| . 
Hence, infinitely many bubbles will be nucleated at dif- 
ferent places and times. The resulting "daughter" bub- 
bles may again contain an asymptotically de Sitter, in- 
finite universe, which again gives rise to infinitely many 
"grand-daughter" bubbles. This picture of eternal self- 
reproduction was called the "recycling universe" [30| . 
Some (or all) of the created bubbles may support a period 
of additional inflation followed by reheating, as shown in 
Fig. El 

In the model of Ref. [30(, there were only two vacua 
which could tunnel into each other. A more recently de- 
veloped paradigm of "string theory landscape" [31[ in- 
volves a very large number of metastable vacua, corre- 
sponding to local minima of an effective potential in field 
space. The value of the potential at each minimum is 
the effective value of the cosmological constant A in the 
corresponding vacuum. Figure 3] shows a phenomenolo- 
gist's view of the "landscape." Vacua with A < do not 
allow any further tunneling 9 and are called "terminal" 
vacua [154| |. while vacua with A > are called "recy- 
clable" since they can tunnel to other vacua with A > 
or A < 0. Bubbles of recyclable vacua will give rise to 
infinitely many nested "daughter" bubbles. A conformal 



6(4>-<f>o). 



(44) 



This is in accord with the fact that Pexit^s) = F^e) 
is a gauge-invariant quantity. Equation (|44p with the 



but 



8 Both processes may be combined in a single scenario [3|_ 
we shall consider them separately for clarity. 

9 Asympto tically flat A = vacua cannot support tunneling Il50l . 
Il5ll . ll52 |: vacu a wit h A < will quickly collapse to a "big crunch" 
singularity [32I . Il53ll . 
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Figure 4: A schematic representation of the "landscape of 
string theory," consisting of a large number of local minima 
of an effective potential. The variable X collectively denotes 
various fields and A is the effective cosmological constant. 
Arrows show possible tunneling transitions between vacua. 




Figure 5: A conformal diagram of the spacetime where self- 
reproduction occurs via bubble nucleation. Regions labeled 
"5" are asymptotically flat (A = 0). 



diagram of the resulting spacetime is outlined in Fig. [3] 
Of course, only a finite number of bubbles can be drawn; 
the bubbles a ctua lly form a fractal structure in a confor- 
mal diagram [l55| . 

A statistical desc riptio n of the "recycling" spacetime 
can be obtained [3(| Il54| by considering a single comov- 
ing worldline x = const that passes through different 
bubbles at different times. (It is implied that the world- 
line is randomly chosen from an ensemble of infinitely 
many such worldlines passing through different points 
x.) Let the index a = 1,...,N label all the available 
types of bubbles. For calculations, it is convenient to 
use the e-folding time r = In a. We are interested in the 
probability f a (r) of passing through a bubble of type a 
at time r. This probability distribution is normalized 
by ^2 a f a = 1; the quantity / Q (r) can be also visual- 
ized as the fraction of the comoving volume occupied by 
bubbles of type a at time r. Denoting by K a p the nucle- 
ation rate for bubbles of type a within bubbles of type /3 
[computed according to Eq. (|46p ]. we write the "master 
equation" describing the evolution of f a (r), 



dr 



^ (- K »0f0 + K /3afa) = ^ M P<xfa, (47) 



where we introduced the auxiliary matrix M a p. Given 
a set of initial conditions / Q (0), one can evolve f a (r) 
according to Eq. (|47| . 

To proceed further, one may now distinguish the fol- 
lowing two cases: Either terminal vacua exist (some j3 



such that K a p = for all a), or all the vacua are re- 
cyclable. (Theory suggests that the former case is more 
probable [53].) If terminal vacua exist, then the late-time 
asymptotic solution can be written as [l54| 



fa(r) « /< 0) + s a e 



(48) 



where fa^ is a constant vector that depends on the initial 
conditions and has nonzero components only in terminal 
vacua, and s a does not depend on initial conditions and 
is an eigenvector of M a p such that J2 a Mp a s a = —qsp, 
q > 0. This solution shows that all comoving volume 
reaches terminal vacua exponentially quickly. (As in the 
case of random-walk inflation, there are infinitely many 
"eternally recycling" points x that never enter any termi- 
nal vacua, but these points form a set of measure zero.) 

If there are no terminal vacua, the solution f a (r) ap- 
proaches a constant distribution [1561 ] . 



lim f a (r) w /, 



(0) 



0. 



(0) 



H a exp 



In this case, the quantities fa^ are independent of initial 
conditions and are interpreted as the fractions of time 
spent by the comoving worldline in bubbles of type a. 

One may adopt another approach and ignore the du- 
ration of time spent by the worldline within each bub- 
ble. Thus, one describes only the sequence of the 
bubbles encount ered along a randomly chosen world- 
line [62|, Il56l . Il57l ] . If the worldline is initially in a bubble 
of type a, then the probability /j,p a of entering the bubble 
of type /3 as the next bubble in the sequence after a is 



7T 

IP 



(49) 



(50) 



M/3o 



E 



(51) 



(For terminal vacua a, we have n ia — and so we may 
define fj,p a — for convenience.) Once again we consider 
landscapes without terminal vacua separately from ter- 
minal landscapes. If there are no terminal vacua, then 
the matrix fi a p is n orma lized, J^p Wo = lj an d is thus 
a stochastic matrix [l58| describing a Markov process of 
choosing the next visited vacuum. The sequence of vis- 
ited vacua is infinite, so one can define the mean fre- 
quency j ( ( ! mcan ) f visiting bubbles of type a. If the prob- 
ability distribution for the first element in the sequence 
is f(Q\ a , then the distribution of vacua after k steps is 
given (in the matrix notation) by the vector 



i(fe) - A* t(o)j 



(52) 



where fi k means the k-th power of the matrix fi = fJ, a p. 
Therefore, the mean frequency of visiting a vacuum a is 
computed as an average of f(k)a over n consecutive steps 
in the limit of large n: 



f< 



n— *oc n 



n 1 n 

n L — ^ v ' n— *oo 77, ' * 



k=l 



(53) 



fc=i 
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(It is proved in the theory of Markov processes that the 
limit j^ mcan ) given by Eq. (|53[) almost surely coincides 

with the mea n frequency of visiting the state a~ see 

e.g. Ref. [l59j |. chapter 5, Theorem 2.1, and Ref. |l60j ] . 
Theorem 3.5.9.) It turns out that the distribution f ( mean ) 
is independent of the initial state f(o) and coincides with 
the distri butio n (I49p found in the continuous-time de- 
scription [156l | . 

If there exist terminal vacua, then almost all sequences 
will have a finite length. The distribution of vacua in a 
randomly chosen sequence is still well-defined and can be 
computed using Eq. (|53|) without the normalizing factor 

n ' 

f (mcan) = (l^rV(o), (54) 

but now the r esul ting distribution depends on the initial 
state f (0 ) [62l.ll57ll. 



III. PREDICTIONS AND MEASURE ISSUES 

As discussed in Sec. II A[ a compelling question in the 
context of eternal inflation is how to make statistical 
predictions of observed parameters. One begins by de- 
termining whether eternal inflation is present in a given 
model. 



A. Presence of eternal inflation 

Since the presence of eternal self-reproduction in mod- 
els of tunneling type is generically certain (unless the 
nucleation rate for terminal vacua is unusually high), in 
this section we restrict our attention to eternal inflation 
of the random- walk type. 

The hallmark of eternal inflation is the unbounded 
growth the total number of independent inflating regions. 
The total proper 3-volume of the inflating domain also 
grows without bound at late times, at least when com- 
puted along hypersurfaces of equal proper time or equal 
scale factor. However, the proper 3-volume is a gauge- 
dependent quantity, and one may construct time gauges 
where the 3-volume decreases wit h time even in an ev- 
erywhere expanding universe [l25l |. The 3-volume of an 
arbitrary family of equal-time hypersurfaces cannot be 
used as a criterion for the presence of eternal inflation. 
However, a weaker criterion is sufficient: Eternal infla- 
tion is present if (and only if) there exists a choice of 
time slicing with an unb ounded growth of the 3-volume 
of inflating domains [l25j . Equivalently, eternal inflation 
is present if a finite comoving volume gives rise to infinite 
physical volume jl38l |. Thus, the presence or absence of 
eternal inflation is a gauge- independent statement. One 
may, of course, use a particular gauge (such as the proper 
time or e-folding time) for calculations, as long as the re- 
sult is known to be gauge- independent. It can be shown 



that eternal inflation is present if and only if the 3-volume 
grows in the e-folding time slicing [l25j . 

The presence of eternal inflation has been analyzed in 
many specific scenarios. For in stance, eternal inflation is 
generic in "chaotic" [H, [U Il6ll ] and "new" 0] inflationary 
models. It is normally sufficient to establish the exis- 
tence of a "diffusion-dominated" regime, that is, a range 
of <j) where the typical amplitude 5<p ~ H of "jumps" is 
larger than the typical change of the field, 0Ai, during 
one Hubble time At = H -1 . For models of scalar-field 
inflation, the condition is 

H 2 > H'. (55) 

Such a range of <f> is present in most slow- roll models of in- 
flation. (For an example of an inflationary scenario whe re 
eternal inflation is generically not present, see Ref. [162| |.) 
A strict formal criterion for the presence of eternal infla- 
tion is the positivity of the largest eigenvalue 7 of the 
operator L + 3H(tfi), as defined in Sec. Ill CI 

The causal structure of the eternally inflating space- 
time and the topology of the reheating surface can be 
visualiz ed u sing the construction of "eternal comoving 
points" [140| | . These are comoving worldlines x = const 
that forever remain within the inflating domain and never 
enter the reheating epoch. These worldlines correspond 
to places where the reheating surface reaches t = 00 in a 
spacetime diagram (see Fig. [2]). It was shown in Ref. [l40j ] 
using topological arguments that the presence of inflat- 
ing domains at arbitrarily late times entails the existence 
of infinitely many such "eternal points." The set of all 
eternal points within a given three-dimensional spacelike 
slice is a measure zero fractal set. The fractal dimension 
of this set can be un derstood as th e fractal dimension of 
the inflating domain [I25lll4fj . ll41| and is invariant under 
any smooth coordinate transformations in the spacetime. 

The existence of eternal points can be used as another 
invariant criterion for the presence of eternal inflation. 
The probability X((f>) of having an eternal point in an 
initial Hubble-size region with field value <fi can be found 
as the solution of a gauge-invariant, nonlinear diffusion 
equation [l40t | 

Dd^X + vd^X - 3H (1 - X) In (1 - X) = 0, (56) 

with zero boundary conditions. Eternal inflation is 
present if there exists a nontrivial solution X(<fi) ^ 
of this equation. 

B. Observer-based measure in eternal inflation 

In theories where observable parameters Xa are dis- 
tributed randomly, one would like to predict the values of 
Xa most likely to be observed by a random (or "typical") 
observer. More generally, one looks for the probability 
distribution P{xa) of observing the values Xa- As dis- 
cussed in Sec. II Al considerations of this type necessarily 
involve some form of the "principle of mediocrity" [l8j . 
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On a more formal level, one needs to construct an en- 
semble of the "possible observers" and to define a proba- 
bility measure on this ensemble. In inflationary cosmol- 
ogy, observers appear only along the reheating surface. If 
eternal inflation is present, the reheating surface contains 
infinitely many causally disconnected and (possibly) sta- 
tistically inequivalent domains. The principal difficulties 
in the probabilistic approach are due to a lack of a natural 
definition of measure on such surfaces. 10 

Existing proposals for an observer-based measure fall 
in two major classes, which may be designated as 
"volume-based" vs. "worldline-based." The difference be- 
tween these classes is in the approach taken to con- 
struct the ensemb le of observers . In the "volume" ap- 
proach [H, EI 51 Ell EH Ell , the ensemble contains 
every observer appearing in the uni verse , at any time o r 
place. In the "wor ldline" approach pU Ell EH Uj5T\ . 
the ensemble consists of observers appearing near a sin- 
gle, randomly selected comoving worldline x = const; 
more generally, an arbitrary timelike geodesic could also 
be used. If the ensemble contains infinitely many ob- 
servers (this is typically the case for volume-based en- 
sembles), a regularization is needed to obtain specific 
probability distributions. Finding and applying suitable 
regularization procedures is a separate technical issue ex- 
plored in Sec. IIII CI and IIII Dl I begin with a general 
discussion of these measure prescriptions. 

A number of previously considered "volume-based" 
measure proposals were f ound to be lacking in one as- 
pect or another [U, EH, EH, EH, EH, Ell, the most 
vexing problem being the dep endence on the choice of 
the time gauge jllll lll2L Il38l |. The requirement of time 
gauge independence is sufficiently important to reject any 
measure proposal that suffers from the gauge ambigu- 
ity. A prescription manifestly free from gauge depen- 
dence is the "spherical cutoff" measure (45|. This pre- 
scription provides unambiguous predictions for models 
of random-walk type eternal inflation if the reheating 
condition = corresponds to a topologically com- 
pact and connected locus in the field space {<f), Xa} (see 
Sec. IIII Cp . For models where the reheating condition is 
met at several disconnected loci in field space (tunncling- 
type eternal inflation belongs to this class), one can 
use the recen tly proposed prescription of "comoving cut- 
off" EM Ess ■ 

Since no other volume-based prescriptions 
are currently considered viable, we refer to the mentioned 
spherical cutoff/comoving cutoff prescriptions simply as 
the "volume-based measure." 

Similarly, existing measure proposals of the "worldline" 
type ap pear t o converge esse ntial l y to a single prescrip- 
tion [H Il57l | (however, see [l7Ct Il7l| for the most re- 
cent developments). We refer to this prescription as the 



"worldline-based measure." 

The main difference between the worldline-based and 
volume-based measures is in their dependence on the ini- 
tial state. When considering the volume-based measure, 
one starts from a finite initial spacelike 3-volume. (Final 
results are insensitive to the choice of this 3-surface in 
spacetime or to its geometry.) The initial state consists 
of the initial values of the fields Xa within the initial vol- 
ume, and possibly a label a corresponding to the type 
of the initial bubble. When considering the worldline- 
based measure, one assumes knowledge of these data at 
the initial point of the worldline. 11 It turns out that the 
volume-based measure always yields results that are in- 
dependent of the initial conditions. This agrees with the 
conce pt o f the "sta tionarity" of a self-reproducing uni- 
verse [ml EH Ell; the universe forgets the initial state 
in the course of eternal self-reproduction. In contrast, 
probabilities obtained using the worldline-based measure 
always depend on the initial state (except for the case of 
a "nonterminal" landscape, i.e. a landscape scenario with- 
out terminal vacua) . A theory of initial conditions is nec- 
essary to obtain a specific prediction from the worldline- 
based measure. 

At this time, there is no consensus as to which of the 
two measures is the physically relevant one. The present 
author is inclined to regard the two measures as reason- 
able answers to two differently posed questions. The first 
question is to determine the probability distribution for 
observed values of Xa, given that the observer is ran- 
domly chosen from all the observers present in the entire 
spacetime. Since we have no knowledge as to the to- 
tal duration of inflation in our past or the total number 
of bubble nucleations preceding the most recent one, it 
appears reasonable to include in the ensemble all the ob- 
servers that will ever appear anywhere in the spacetime. 
The answer to the first question is thus provided by the 
volume-based measure. 

The second question is posed in a rather different man- 
ner. In the context of tunneling-type eternal inflation, 
upon discovering the type of our bubble we may wish to 
leave a message to a future civilization that may arise in 
our future after an unspecified number of nested bubble 
nucleations. The analogous situation in the context of 
random-walk type inflation is a hypothetical observer lo- 
cated within an inflating region of spacetime who wishes 
to communicate with future civilizations that will even- 
tually appear when inflation is over. The only available 
means of communication is leaving information on paper 
in a sealed box. The message might contain the probabil- 
ity distribution P(xa) for parameters Xa that we expect 
the future civilization to observe. In this case, the initial 



To avoid confusion, let us note that the recent work Il6cfl pro- 
poses a measure in the phase space of trajectories rather than an 
observer-based measure in the sense discussed here. 



1 Naturally, it is assumed that the initial state is in the self- 
reproduction regime: For random-walk type models, the initial 
3-volume is undergoing inflation rather than reheating and, for 
tunneling models, the initial 3-volume is not situated within a 
terminal bubble. 
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Figure 6: Illustrative inflationary potential with a flat self- 
reproduction regime (ft\ < (ft < <j>2 and deterministic regimes 
<j>i 1} < (ft < 0i and (j> 2 < (ft< 4 2) . 

state is known at the time of writing the message. It is 
clear that the box can be discovered only by future ob- 
servers near its worldline. It is then natural to choose 
the ensemble of observers appearing along this worldline. 
Starting from the known initial state, one would then 
compute P(Xa) according to the worldline-based mea- 
sure. 

Calculations using the worldline-based measure usu- 
ally do not require regularization (except for the case of 
nonterminal landscape) because the worldline-based en- 
semble of observers is almost surely finite [157j | . For in- 
stance, in random-walk type models the worldline-based 
measure predicts the exit probability distribution p e xitj 
which can be computed by solving a suitable differential 
equation [see Eq. (|44|)], However, the ensemble used in 
the volume-based measure is infinite and requires a reg- 
ularization. Known regularization methods are reviewed 
in Sections UTTC] and niDl 

A simple toy model [18( where the predictions of the 
volume-based measure can be obtained analytically is a 
slow-roll scenario with a potential shown in Fig. [6] The 
potential is flat in the range (ft\ < (ft < 4>2 where the 
evolution is diffusion-dominated, while the evolution of 
regions with (ft > (ft>2 or (ft < (fti is completely determin- 
istic (fluctuation- free). It is assumed that the diffusion- 
dominated range (fti < (ft < <\>i is sufficiently wide to cause 
eternal self- reproduction of inflating regions. There are 
two thermalization points, (ft — (fti 1 ^ and (ft = (ft* , which 
may be associated to different types of true vacuum and 
thus to different observed values of cosmological param- 
eters. The question is to compare the volumes V\ and V2 
of regions thermalized into these two vacua. Since there 
is an infinite volume thermalized into either vacuum, one 
looks for the volume ratio Vi/V2- 

The potential is symmetric in the range (ft\ < (ft < <p2 1 
so it is natural to assume that Hubble-size regions exiting 
the self-reproduction regime at (ft — (fti and at (ft = (f>2 are 
equally abundant. Since the evolution within the ranges 

(fti 1 ^ < (ft < (fti and 4>2 < (ft < (ft^ is deterministic, the 
regions exiting the self-reproduction regime at (ft — (fti or 
(ft = (ft>2 will be expanded by fixed amounts of e-foldings, 
which we may denote Ni and N2 respectively, 

Therefore the volume of regions thermalized at (ft — (fti^ , 




-V 



Figure 7: A 1+1-dimensional slice of spacetime in a two-field 
inflationary model (numerical simulation in Ref. [3). Shades 
denote different values of the field x> which takes values in 
the periodically identified interval [0, 2n//3]. The white re- 
gion represents the thermalized domain. The boundary of 
the thermalized domain is the reheating surface (cf. Fig. [2]| , 
which contains all the possible values of the field \. 



where j = 1, 2, will be increased by the factors exp(3AT,). 
Hence, the volume ratio is 

£ = exp (37V! - 3N 2 ) . (58) 

V2 



C. Regularization for a single reheating surface 

The task at hand is to define a measure that ascribes 
equal weight to each observer ever appearing anywhere 
in the universe. As discussed in Sec. II A[ it is sufficient 
to construct a measure V(xa) of the 3-volume along the 
reheating surface. The volume-based measure P{xa) will 
then be given by Eq. (f5]). In the presence of eternal 
inflation, the proper 3-volume of the reheating surface 
diverges even when we limit the spacetime domain un- 
der consideration to the comoving future of a finite ini- 
tial spacelike 3-volume. Therefore, the reheating surface 
needs to be regularized. 

In this section we consider the case when the reheat- 
ing condition is met at a topologically compact and con- 
nected locus in the configuration space {(ft,Xa\- In this 
case, every connected component of the reheating sur- 
face in spacetime will contain all the possible values of 
the fields Xai an d all such connected pieces are statisti- 
cally equivalent. Hence, it suffices to consider a single 
connected piece of the reheating surface. A situation of 
this type is illustrated in Fig. [71 A sketch of the random 
walk in configuration space is shown in Fig. [5] 

A simple regularization scheme is known as the "equal- 
time cutoff." One considers the part of the reheating sur- 
face formed before a fixed time i max ; that part is finite 
as long as i max is finite. Then one can compute the dis- 
tribution of the quantities of interest within that part of 
the reheating surface. Subsequently, one takes the limit 
t max — > 00. The resulting distribution can be found from 



14 




Figure 8: A random walk in configuration space for a two- 
field inflationary model considered in Ref. [19( ]. The center of 
the field space is a diffusion-dominated regime. The reheating 
condition, (f> — </)*, selects a compact and connected region (a 
circle) in configuration space. The problem is to determine 
the volume- weighted probability distribution for the values of 
X at reheating. 

the solution Py ((f)) of the "stationary" FP equation, 

[L + 3H(<f>)]P™=>yP<?\ (59) 

with the largest eigenvalue 7 (see Sec. Ill CI) . However, 
both the eigenvalue 7 and the distribution Py\tfi) de- 
pend rather sensitively on the choice of the equal-time 
hypersurfaces. Since there appears to be no preferred 
choice of the cutoff hypersurfaces in spacetime, the equal- 
time cutoff cannot serv e as an unbiased measure. Also, it 
was shown in Ref. [1251 ] that the unbiased result (|58|) can- 
not be obtained via an equal-time cutoff with any choice 
of the time gauge. 

The "spherical cutoff" measure prescription [45[ reg- 
ularizes the reheating surface in a different way. A fi- 
nite region within the reheating surface is selected as 
a spherical region of radius R around a randomly cho- 
sen center point. (Since the reheating 3-surface is space- 
like, the distance between points can be calculated as 
the length of the shortest path within the reheating 3- 
surface.) Then the distribution of the quantities of in- 
terest is computed within the spherical region. Subse- 
quently, the limit R — > 00 is evaluated. Since every por- 
tion of the reheating surface is statistically the same, the 
results are independent of the choice of the center point. 
The spherical cutoff is gauge-invariant since it is formu- 
lated entirely in terms of the intrinsic properties of the 
reheating surface. 

While the spherical cutoff prescription successfully 
solves the problem of regularization, there is no univer- 
sally applicable analytic formula for the resulting distri- 
bution. Application of the spherical cutoff to general 
models of random-walk inflation requires a direct numer- 
ical simulation of the stochastic field dynamics in the 
inflationary spacetime. Su ch simulations were reported 
in Refs. [H, Ill2l Il4ll Il72l | and used the Langevin equa- 
tion ([9]) with a specific stochastic ansatz for the noise 
field N(t,x). 

Apart from numerical simulations, the results of the 
spherical cutoff method may be obtained analytically in 



a certain class of models . One such case is a multi- 
field model where the potential V(4>, Xa) is independent 
of Xa within the range of <f> where the "diffusion" in <p 
dominates. Then the distribution in Xa is flat when field 
</> exits the regime of self-reproduction and resumes the 
deterministic slow-roll evolution. One can derive a gauge- 
invariant Fokker-Planck equation for the volume distri- 
bution Pv(4>, Xa)i using <j) as the time variable [19j , 

d*Pv = dxx (-Pv) ~ d x (^Pv) + —Pv, (60) 
W / \ v 4> J v 4> 

where D, v^, and v x are the coefficients of the FP equa- 
tion pT|) . This equation is valid for the range of <f> where 
the evolution of <j> is free of fluctuations. By solving 
Eq. (|5D|) . one can calculate the volume-based distribu- 
tion of Xa predicted by the spherical cutoff method as 
Pv (4>*, x)- Note that the mentioned restriction on the 
potential V(</>, Xa) is important. In general, the field 4> 
cannot be used as the time variable since the surfaces 
of constant <j> are not everywhere spacelike due to large 
fluctuations of <fi in the diffusion-dominated regime. 

D. Regularization for multiple types of reheating 
surfaces 

Let us begin by considering a simpler example: an in- 
flationary scenario with an asymmetric slow-roll poten- 
tial V(4>) having two minima 4>^\ 4& ■ This scenario 
has two possibilities for thermalization, possibly differ- 
ing in the observable parameters Xa- More generally, one 
may consider a scenario with n different minima of the 

potential, possibly representing n distinct reheating sce- 

(i) 

narios. It is important that the minima <pl , j = 1, ...,n 
are topologically disconnected in the configuration space. 
This precludes the possibility that different minima are 
reached within one connected component of the reheat- 
ing hypersurface in spacetime. Additionally, the fields Xa 
may fluctuate across each connected component in a way 
that depends on the minimum j. Thus, the distribution 
P(Xa]j) of the fields Xa at each connected component of 
the reheating surface may depend on j. In other words, 
the different components of the reheating surface may 
be statistically inequivalent with respect to the distribu- 
tion of Xa ° n them. To use the volume-based measure 
for making predictions in such models, one needs a regu- 
larization method that is applicable to situations with a 
large number of disconnected and statistically inequiva- 
lent components of the reheating hypersurface. 

In such situations, the spherical cutoff prescription (see 
Sec. IIII C|) yields only the distribution of Xa across one 
connected component, since the sphere of a finite radius 
R will never reach any other components of the reheat- 
ing surface. Therefore, the spherical cutoff needs to be 
supplemented by a "weighting" prescription, which would 
assign a weights pj to the minimum labeled j. In sce- 
narios of tunneling type, such as the string theory land- 
scape, observers may find themselves in bubbles of types 
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a = 1,...,N. Again, a weighting prescription is needed 
to determine the probabilities p a of being in a bubble of 
type a. 

Two different weighting prescript ions have been for- 
mulated, using the volume-based [l5l [HI and the 
worldline-based approach |l57| . respectively. The first 
prescription is called the "comoving cutoff" while the sec- 
ond the "worldline" or the "holographic" cutoff. For clar- 
ity, we illustrate these weighting prescriptions on mod- 
els of tunneling type where infinitely many nested bub- 
bles of types a — 1,...N are created and where some 
bubbles are "terminal," i.e. contain no "daughter" bub- 
bles. In the volume-based approach, each bubble receives 
equal weight in the ensemble; in the worldline-based ap- 
proach, only bubbles intersected by a selected worldline 
are counted and given equal weight. Let us now examine 
these two prescriptions in more detail. 

Since the set of all bubbles is infinite, one needs to 
perform a "regularization," that is, one needs to select a 
very large but finite subset of bubbles. The weight p a 
will be calculated as the fraction of bubbles of type a 
within the selected subset; then the number of bubbles 
in the subset will be taken to infinity. Technically, the 
two prescriptions differ in the details of the regulariza- 
tion. The difference between the two prescriptions can 
be understood pictorially (Fig. [5]). In a spacetime dia- 
gram, one draws a finite number of timelike comoving 
geodesic worldlines emitted from an initial 3-surface to- 
wards the future. (It can be shown that the results are 
independent of the choice of these lines, as long as that 
choi ce is uncorrelated with the bubble nucleation pro- 
cess [l64j].) Each of these lines will intersect only a finite 
number of bubbles, since the final state of any worldline is 
(almost surely) a terminal bubble. The subset of bubbles 
needed for the regularization procedure is defined as the 
set of all bubbles intersected by at least one line. At this 
point, the volume-based approach assigns equal weight 
to each bubble in the subset, while the worldline-based 
approach assigns equal weight to each bubble along each 
worldline. As a result, the volume-based measure counts 
each bubble in the subset only once, while the worldline- 
based measure counts each bubble as many times as it 
is intersected by some worldlines. After determining the 
weights p a by counting the bubbles as described, one in- 
creases the number of worldlines to infinity and evaluates 
the limit values of p a . 

It is clear that the volume-based measure represents 
the counting of bubbles in the entire universe, and it is 
then appropriate that each bubble is being counted only 
once. On the other hand, the worldline-based measure 
counts bubbles occurring along a single worldline, ignor- 
ing the bubbles produced in other parts of the universe 
and introducing an unavoidable bias due to the initial 
conditions at the starting point of the worldline. 



Explicit formulas for p a were derived for a tunneling- 
type scenario (with terminal vacua) in the volume-based 




Figure 9: Weighting prescriptions for models of tunneling 
type. Shaded regions are bubbles of different types. Dashed 
vertical lines represent randomly chosen comoving geodesies 
used to define a finite subset of bubbles. 



approach [15 



(61) 



where K a p is the matrix of nucleation rates, q and s a 
are the quantities defined by Eq. and Hp is the 

Hubble parameter in bubbles of type (3. The expressions 
for p a obtained from the worldline-based measure are 
given by Eq. (f5"4")l . As we have noted before, the volume- 
based measure assigns weights p a that are independent 
of initial conditions, while the weights obtained from the 
worldline-based measure depend sensitively on the type 
of bubble where the counting begins. 

In the case of a non-terminal landscape, both the 
volume-based and the world-line based measures give 
identical results for p a , which coincide with the mea n 
frequency (T3!)j) of visiting a bubble of type a (l56l 1 157j ] . 

With the weighting prescriptions just described, the 
volume-based and the worldline-based measure propos- 
als can be considered complete. In other words, we have 
two alternative prescriptions that can be applied (in prin- 
ciple) to arbitrary models of random-walk or tunneling- 
type eternal inflation. Further research is needed to reach 
a definite conclusion concerning the viability of these 
measure prescriptions. 
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